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Abstract 



1 Introduction 



In matrix model proposals for type IIB superstring theory and M-theory CD HI, describing gravity 
by matrices is one of the important subjects. Gravity is not prepared at formulations of these matrix 
models, but it must be encoded in their matrix degrees of freedom. Several years ago, a new inter- 
pretation of large N matrix was proposed by Hanada, Kawai and Kimura (3j |4| ; large N matrices 
can describe covariant derivatives in curved space. They applied this new interpretation (we call 
this interpretation "HKK interpretation" or simply "HKK" for later convenience) to the IKKT matrix 
model [1], and derived the Einstein equation in empty space from the classical matrix equation of 
motion. Several studies have been done to explore this attractive proposal f5l 161171 151 |9]R 

In the HKK interpretation, large ,/V matrices contain not only vielbein and spin connection but 
also infinitely many other degrees of freedom with various higher spins. If one consider gravity de- 
rived by HKK without these higher spin fields, this is regarded as an effective theory in nearly empty 
space. Even under such simplification, however, this effective theory has more degrees of freedom 
than Einstein gravity. In HKK, the vielbein and the spin connection encoded in large ,/V matrices are 
independent variables, then the matrix description of gravity gives rise to torsion degrees of freedom 
unless we impose some torsion free conditions by hand. Therefore IKKT matrix model provides a 
torsion gravity. This arouses our curiosity on the role of torsion in the gravity from the matrix model. 
The classical field equation would be a good way to capture its semi-classical dynamics. Gravity 
equations with torsion from several bosonic IKKT-type matrix models have been written down in 
(51 . Our purpose in this paper is to investigate classical solutions of the gravity equation with torsion, 
from the bosonic IKKT model, under simple settings of torsion. 

There are two remarks: 

• A generalization of Einstein gravity to including torsion has been known as the Poincare gauge 
theory of gravity (PGT), and has been studied extensively (see ifTTl [121 for reviews). It is 
described by a gravity Lagrangian which includes bi-linear terms of torsion. Vierbein and spin 
connection in this Lagrangian are independent variables. On the other hand, in this paper we 
will study a gravity equation with torsion which is derived from the HKK interpretation of 
a matrix model equation of motion, which is different from what is derived from the PGT. 
Although the matrix equation of motion is derived from a matrix model action, the gravity 
equation of HKK seems difficult to derive by a variation of any gravitational Lagrangian. 

• In (5J, parts of torsion degrees of freedom were identified with scalar and 2nd rank antisym- 
metric tensor fields, which are analogous to the massless fields of string theory. In this paper 
we do not always follow such an interpretation but rather keep broader possibility for solutions. 

3 Other similar attempts based on non-commutativity can be found in 1101 
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The organization of this paper is as follows: In section 2, we briefly describe the gravity equation 
with torsion from the HKK interpretation and fix our notation. The original proposal was done for 
a matrix model with an Euclidean signature. After writing down the field equation with Euclidean 
signature, we rotate it to one with a Lorentzian signature. In section 3, we consider time dependent 
solutions under homogeneous-isotropic ansatz. In section 4, static and spherical symmetric solutions 
are studied. In addition to the analytic method, we examine numerical computation to find solutions. 
Finally we summarize and discuss our results in section 5. 

2 Gravity Equation from Large N Matrix Model 

Large N matrices can describe covariant derivatives in curved spaces (3HH. In this interpretation, the 
large N matrices are regarded as a map between sections of a fiber bundle over a curved manifold 
whose fiber is the vector space of structure group Spin(J). Formally differential operators may be 
expanded by a power series of covariant derivatives and Lorentz generators. Associated with this 
expansion, various fields: U(l) gauge field, vielbein, spin connection and other infinitely many 
higher spin fields appear. Here we concentrate ourselves on the degrees of freedom associated with 
the vielbein and spin connection. In this situation the HKK interpretation is simply stated as 

A {a )=R { a) a (g- i )V a , (2-1) 

where index (a) = 1 , 2, ...,d labels N xN Hermitian matrices, on the other hand, another index a is the 
local Lorentz indices in J-dimensional curved manifold. R^afis 1 ) i s m the vector representation 
of the local Lorentz group Spin(ci) whose elements are denoted by g. We describe curved space 
covariant derivative V b by the vielbein e M a and the spin connection iO a bc as 

y a = e f ' a d IJ + (O a bc O bc , (2.2) 

where /u is the curved space index and O bc is the Lorentz generator. 
Commutator of these covariant derivatives gives 

[V a ,V*] = -T c , th W c + R ab cd o cd . (2.3) 

Here T c xlb and R a b cd are defended as 

-T\ cl b = ep e @ a J'b-dbJ , a)+(Bab e -aiba e , 
R ab cd = d a (O b cd -d b (O c cd + 03 a c e (O b d e -(O b c e 03 a d e (2.4) 

where d a = e^ad^. We do not impose torsion free condition by hand. Thus the torsion T c ab is 
non-vanishing and e^a and (£> a bc are independent of each other. 



2 



Next, to obtain a gravity equation, we apply this interpretation to the IKKT-type matrix model 
(d = 10 conesponds to the IKKT model). The bosonic part of the matrix model action is 

- iTr[A (a) ,A (fo) ][A (c) ,A w ]5W( c )5WW. (2.5) 
Matrix equation of motion from this action becomes 

[A^,[A {a) ,A {b) }}=0. (2.6) 

We now apply the HKK interpretation (12.11 ) to (12.61) . With a help of a property (68) in 0, this 
equation becomes 

[V u ,[V a ,V,]]=0. (2.7) 

From this equation we obtain following equations 

Rab + V c T hiCa + T™ a T h , pil = 0, 

V a Rabcd + T™ b R pq cd = 0, (2.8) 

where R ab = R aeb e . 

From a matrix Jacobi identity 

[v w , [V , V (c) ]] + [V {b) , [V (c) , v (fl) ]] + [v (c) , [V (a) , V (b) ]} = 0, 

we have another set of equations: 

- ^( a T d .bc) ~ R (bca) d + T e (,bcT d ,a)e = 0, 

V(aRbc) d e -T f Aab R c)f d e = 0. (2.9) 

These equations are identical with Bianchi identities of T c ^ ab and R ab cd . Thus (12.91 ) is automatically 
satisfied according to the definitions (I2.4I ). See the Appendix as a proof. 

Equations (12.81 ) are derived from IKKT-type matrix model with Euclidean signature. Having a 
gravity equation once, we rotate it to have Lorentzian signature in order to discuss classical solutions 
which describe manifolds with Lorentzian signature. Then, (12.81 ) with Minkowski signature are basic 
equations which we will use from the next section. 

3 Time Dependent Solutions with Homogeneity and Isotropy 

In this section we consider time dependent solutions with assumption of homogeneity and isotropy. 
General ansatz for the spin connection under this assumption was already proposed in PGT |[T3l . 
Here we do not use this most general ansatz. Rather we shall consider simple settings as possible, 
and examine the role of the torsion in our classical gravity. 
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3.1 Time dependent solutions with single scale factor 

We shall consider (n + 1) -dimensional homogeneous and isotropic spacetime. For the metric and the 
spin connection we adapt following ansatz: 



ds 1 = -dt 2 + a(t) 2 ^ dxf, 

(Oi 1 = (02 2 o = ••• = ® n "o = X(t), (others) =0. (3.1) 
Under this ansatz non-vanishing components of torsion are 



(3.2) 



Substituting these forms into (12.8I ). we obtain following equations: 



X+X 2 = 0, (3.3) 
2X-H + (2n-\)X 2 -nHX + H 2 = 0, (3.4) 
(d t + nX-H)(X + HX) - (n - l)X 3 = 0, (3.5) 



where (13.31 ) and (13.41 ) come from the first equation of (12.8I ). and (13.51 ) comes from the second equation 
of (12.81) . We have three equations here, but it turns out that (13.51 ) can be written using (13.31 ) and (I3.4I ). 
Therefore we consider (13.31 ) and (13.4b as independent equations to determine X and H. 

First, we find a solution of (I3.3I )- (I3.5I ) with following form: 

X = 0, // = —!—, (3.6) 
a-t 

with an integral constant a. The scale factor becomes a(t) °< (a — t)~ l , so it blows up as time with 
positive acceleration i.e. a > 0. In this solution all quantities calculated by Riemann curvatures 
vanish because X is zero. Expansion of the space by a(t) is purely supported by the torsion. 

Next we consider X / solution. We can solve (13.31 ) by X = (f3 + t)~ l where [3 is an integral 
constant. Substitute it into (13.41 ). we have 



r2 n_ _2n-3 



H = H 1 - ir —H+ - : (3.7) 



A solution of (13.71 ) can be written by 



F-jlj, (3.8) 



and m is determined by second order equation 

u 2 - (n- 1)m + (2n-3) =0. (3.9) 
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There are real solutions of u if n = 1 or n > 9. The scale factor becomes a{t) <=c (t + (3)" and a > is 
possible if n > 9. 

In the case of 1 < n < 9, there is no simple analytic solution. However, ( 13.71 ) tells us that its 
solution approaches to (13.61 ) as t increases. Thus finally it blows up with positive acceleration. 



Every three cases which we have considered implys that introducing torsion gives expansion with 
positive acceleration to the universe. 

3.2 Time dependent solutions with two scale factors 

We shall consider (n + m+ 1) -dimensional homogeneous and isotropic spacetime with two different 
scale factors. For metric and spin connections we adapt the following ansatz: 

n n+m 

ds 2 = -dt 2 + a(t) 2 Y,dxf + b(t) 2 £ dxj, 

1=1 I=n+l 

ra/o = »2 2 o = • • • = co/o = *(0> co„ +1 ' !+1 = aW l+2 o = - = <W +OT o = Y(t), 
(others) = 0. (3.10) 

Under this ansatz non-vanishing components of torsion are 

T\ i0 =X-H 1 , T I J0 = Y-H 2 , (h y = ^, H 2 = ^\. (3.11) 



By substituting these forms into (12.81 ). we obtain following equations: 

n{X+X 2 )+m(Y + Y 2 ) = 0, (3.12) 

2X - H x + (2n - 1 )X 2 + 2mXF - (nX + mF)ffi + H\ = 0, (3. 13) 

2Y - H 2 + (2m - 1 )Y 2 + 2nXY - (nX + mY )H 2 +H% = 0, (3. 14) 

X+XHx + (nX + mY)X + (nX + mY - H x )XH X - (n- l)X 3 - mX 2 Y = 0, (3.15) 

Y + YH 2 + (nX + mY)Y + (nZ + mY - H 2 )YH 2 - (m - 1)F 3 - nXY 2 = 0, (3.16) 

where (13.121 ), ( 13.131 ) and ( 13.141 ) come from the first equation of (I2.8I ). On the other hand, ( 13.151 ) and 
(13.161 ) come from the second equation of (12.81) . We have five equations here, but it turns out that 
one of these five is not an independent equation^. Therefore we regard (13.12l )- (13.15l) as independent 
equations to determine X,Y,H\ and H 2 . In the paper (51, torsion components T'jo and T 1 jq are 
identified with a time derivative of some scalar field. Thus if we follow this proposal, additional 
constraint X — Hi =Y — H 2 must be imposed. But here we do not impose this constraint and keep a 
broader possibility for the solution. 



4 By using |gH=0, one can see nX x $3A3j +mY x (|3l4t =nx ( IXTsT l +m x 1X161 ). 
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The same as in the previous subsection, the field equations above have the following curvature- 
less solution: 

X = Y = 0, Hi = ^—,H 2 = ^—, (3.17) 

where a\ 2 are integral constants. 

It is difficult to have an analytic result in the case of X,Y ^ 0. Here we consider linear approxi- 
mation of the equations d3.12| )-( [3~T5T ): 

X=0, nX + mY = 0, 
2X-H 1= 0, 2Y-H 2 = 0. (3.18) 

These are easily solved by 

X = pf + P', Y = p? + p", H x = 2p? + y H 2 = p?+y, (3.19) 

m m 

with integral constants P, P', P" and y y Let y = y = 0, for example. Then we have a{t) exp(pV 2 ) 
and b{t) oc exp(— ^t 2 ). If a{t) gives expansion or shrinking in*' directions, then b{t) gives shrinking 
or expansion in y 1 directions. Acceleration a is positive if (3 > 0. Acceleration b is negative at early 
time and becomes positive after t > m/2np\ Linear approximation may be valid during t < 1/20. 
Thus b can be positive if m/n < 1 , within the linear approximation. 

For large X,Y,Hi and H 2 we have to take care of full non-linearity of the equations (I3.12I )- (I3.15I >. 
Numerical computation is effective to do it. Depending on initial conditions for (X ,Y,H\,H 2 ), we 
may experience various periods of evolution; a(t),b(t) give expansion or shrinking with or without 
acceleration. We left such detailed numerical analysis to future work. 

4 Static Solutions with Spherical Symmetry 

In this section we study solutions which describe static solutions with spherical symmetry. We re- 
strict our interest to (3+l)-dimensional spacetime. The same as in the previous sections, we do not 
intend to exhaust the most general ansatz which respects spherical symmetry for the spin connection. 
Instead, we set a simple ansatz in order to investigate the role of torsion. 

We shall adopt the following ansatz for the metric and the spin connection: 

ds 2 = -F{r) 2 dt 2 + G(r) 2 dr 2 + r 2 (dQ 2 + sin 2 6^ 2 ), 

COS 

C0oi= A(r), (0 2 2 i =(0 3 3 i =B(r), (0 3 3 2 = ^-, (others) =0. (4.1) 

rsmB 

Under this ansatz the non-vanishing components of the torsion are 

r° 01 =A-D, T 2 . 2l =T\^=B-E, (d = ^,E = ^]. (4.2) 
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Substituting these forms into ( 12.81 ), we obtain the following equations: 



— (nA' + mB') + nA 2 + mB 2 = 0, 
G 



(2A' - D') +AD + mAB + mB{A -D) + {A- D) 2 = 0, 
(2B' -E')+BE + nAB + nA(B-E) + (B-E) 2 + P = 0, 



d,+A + 2B-D 
d,+A + 2B-E 



A'+AD 

B' + BE 
1 



mAB = 0, 



nA 2 B + Q = 0, 



1 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 



We introduced n = \,m = 2 to write down these equations in more symmetric form. Equations (14.3b . 
( 14.41 ) and (14.51 ) come from the first equation of ( I2.8I ). On the other hand, ( 14.61 ) and (14.71 ) come from 



the second equation of 
an independent equatior 



( 2.8b . There are five equations, however, it turns out that one of them is not 
Thus we use four equations d4.3l )-( |4~6l ) to determine A, B,D, and E. In the 
paper (H, r°,oi, j ,21, and r 3 ,3i are identified with a derivative of some scalar field. If we follow 
their proposal then additional constraint A — D = B — E must be imposed. Again we do not impose 
this constraint here in order to keep a broader possibility for the solutions. 



4.1 Analytic approach 

In the torsion-less case, the Schwarzschild spacetime F 2 = \/G 2 = 1 — co/r is a solution. It corre- 
sponds to 



A=D 



= , B = E= l -Jl- C ^. 



(4.8) 



Similar to the previous sections, we can find a curvature-less solution A = B = 0. The equations 
for D and E are 

^D'-D 2 = 0, ^E'-E 2 + \=0. (4.9) 
G G r L 

Recall that D = F' /FG and E = 1/rG, then we find two different solutions of (I4.9I ). The first solution 

is 



with integral constants c,c' , and c"Q The second solution is 



e' 

F = , G = —=, (A=B = 0) 

e-r ^2 



(4.10) 



(4.11) 



5 One can see n x [|4~6ll +m x fi~7l) = 0, using B3t . fi~4t . i POt and ^3,-P3t. 
6 1/F can be written as an elliptic integral. 
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with integral constants £ and e'. These are not an asymptotically flat solution. 

It is difficult to have analytic solution for A, B ^ with torsion. Here we consider a linear 
approximation of the equations d4.3l )-( |4~6l ): 



LdJ±d rA }=0, L(A' + 2B')=0, 



-(24'-D')=0, ~(2B' -E')=0. (4.12) 
G G 

These equations may be valid for the large r region. From these equations, we have 

—A' = —8, -B' = 28, -D' = -25, -£' = 45. (§ = const.) (4.13) 
G G G G 

Recalling E = 1/ rG, then the solution is 

A = -8 J' G{s)ds, B = 28 J" G(s)ds, 



lo 



gF = -8^J r G(s)ds^ +8', J-_ = (88logr + 8"). (4.14) 

Again, this does not give asymptotically flat spacetime. E behaves as log r, so this linear approxima- 
tion breaks down not only in the small r region but also in a very large r region. 

An asymptotically flat solution can be constructed by a formal power series expansion. To do 
this we rewrite (14.3I )- (14.6I) as first order differential equations: 



dA 


C 




dr 


Tr' 




dB 


\C+A 2 + 2B 2 




~dr~ 


2 Er 




dC 


- (3A + 2B)C - 4A 2 B + 2AB 2 


-A 3 


dr 


Er 




dD 


2C+AD + 2B(A -D) + (A- 


Df 


dr 


Er 




dE 


-C-A 2 + 2BE + (B-E) 2 - 


l/r 2 


dr 


Er 





(4.15) 

Coefficients of formal l/r expansion can be determined by (14- 15b - The result is 

a a(3a — 2b) , A . „ d a 2 + 3ad — 2bd + d 2 . 4 . 
A = 7 2+ r 3 ' +0{r-% D=- 2+ -3 + 0{r-% 

„ 1 b a 2 -Sab + 5b 2 1 . 4 , la 13a 2 + I2ab + b 2 1 . 4n 

B = - r ~7 2+ 4 7 + ° {r } ' E = -r~7 2 4 7 + ° {r ] - 

(4.16) 

Constants a,b,d parametrize the solution. From this expression we obtain the metric 

F - 1 _j_^+4a- M - M ) +o(r _ 3)j 

1 a 13a 2 - \2ab + b 2 1 , 

- = 1 2+0(r ~ 3 ), (4.17) 

& r 4 r z 
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and the torsion 



a — d 2a 2 — 2ab — 3ad + 2bd + d 2 , dx 
A-D = — + 3 +Or" 4 



B-E = «=± + la2 - l °? + 3b2 ^ + 0(r-\ (4.18) 

Although we can formally continue to construct the power series solution to higher order, it is not 
clear whether the event horizon does appear or not. From the next subsection we examine a nu- 
merical integration to treat full non-linearity of equations and discuss the existence of horizon. In 
the numerical approach, a generic boundary condition does not give an asymptotically flat solution. 
Probably the formal power series (14.161 ) is not a convergent series and cannot be a full solution. We 
may use it as an approximately solution which describes the asymptotic form of a spacetime with 
some matter in the inner region. 

4.2 Numerical approach 

A numerical approach is effective to treat full non-linearity of coupled differential equations. For 
convenience of the analysis, we change variables in differential equations (14.151 ) as 

1 

r= -, 
u 

A^-u 2 A, Dh^u 2 D, C*-^u 3 C, B^-uB, E^-uE. 

After this, changing (14.151 ) becomes 

dA 1 / C 

— = — 2A + - 
du u V t 



dB 1 / B 2 \ \C+A 2 u 



du u \ E ) 2 E 

dC 1 ( 2BC-2AB 2 \ 3AC + 4A 2 B+A 3 u 
-3A H + ■ 



du u \ E J E 

dD 1/ 2C + AAB-2BD\ AD+lA-D) 2 
— 2D+ + 



du u \ E J E 

dE If (B-E) 2 -l\ C-2AB+AE+A 2 u 

— = -lE + i ^ + . (4.19) 

du u \ E J E 

We used the "desolve" function of Maple to obtain a numerical solution, imposing a boundary condi- 
tion at u = -7-U. We have found four different classes of the numerical solutions which are character- 

n 

ized by singularitieo Some of them seem to correspond to the event horizon. We adopt an intuitive 
criterion for the horizon; if numerically F(r)/G(r) = 0, then it is identified with the existence of 
the horizon, because F{r)/G(r) is the speed of light seen by an observer at infinity in the case of 
Schwarzschild spacetime. Below we list up solutions we found: 



7 Here what singularities mean should be understood in a numerical sense, though we expect that it relates to a non- 
analytic property of solution. 
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• Type I : There are two singularities at u = and w = w/,. The solution exists between < u < w/,. 
The point u = u/, is identified with the event horizon. Around u = 0, the solution blows up with 
power laws. 

• Type II : There are two singularities which are identified with horizons. The solution exists 
between these two horizons. 

• Type III: There is a singularity at u = which is identified with the horizon. The solution 
exists in u/, < u. 

• Type IV : There is a singularity at u = 0, but this is not a horizon, The solution exists in < u. 
Around u = 0, it blows up with power laws. 

The curvature tensor R a b cd and the scalar curvature R = R a b ah can diverge at these horizons. However, 
it could be canceled by divergence of the torsion. According to the field equation (12.81 ), a particular 
scalar combination: R + V c T c \ ca + T p,qa T p ^ qa is always zero except at positions of sources. In this 
sense the divergences of the curvature tensor and the torsion at these horizons do not immediately 
mean the singularity of the solution. 

In next subsection, we discuss these solutions in more detail. 

4.3 Details of numerical solutions 
type I solution 

An example of type I solution is given by the boundary condition: 



Figure Q] shows that there are two singularities for solutions. One is located at u < 10~ 20 . We think 
that the real position of this singularity is u = 0, and it is caused from power law behavior of the 
solutions around u ~ 0. Another singularity is located at = 21.4. Figure [2] shows that near u = m/,, 
D(u) and E(u) behave like 




(4.20) 



D(u) ~ D h (u h - u) p , p<0, 
E(u) ~ Eh(uh — u) q , < q < 1 



(4.21) 



with some positive number D/, and £"/, 



Recalling D(r) = and E{r) = ^, then 




(4.22) 



X 



Note if y = (x - x) q , then y" < for < q < 1. 
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Figure 1: Type I solution: Log-Log plot of D and E. Horizontal line is x = logw. 



0.1- 



0.0B 



0.06 



0.04 



0.02- 




ure 2: Type I solution: plots of D and E near the singularity U] x = 21.4. Horizontal lines are 
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Figure 3: Type II solution: plots of D and E. Horizontal lines are u. 



In the case of \p + q\ — 1 ^ 0, where C is some positive numerical factor. If \p + q\ — 1 > 0, then 
F — > as u — > Uh, while F becomes finite at u = Uh is \p + q\ — 1 < 0. Then F/G = at w/,. In the 
case of + = 1, we have F ~ log(«/, — u). In this case we evaluate F/G as 



F 
G 



U — >U)x 



G, 



U—fUfa 



oc (U h -U) q \ u ^ llh =0. 



u— mi, 



We used < q < 1 in the last step. Now F/G = at w = Uf, has been established. Then we conclude 
Uh is the event horizon according to our criterion. The type I solution describes the outside of a black 
hole in an asymptotically non-flat spacetime. 



type II solution 

An example of the type II solution is given by boundary condition: 

x / 1 —2 1 1 \ , 

(A,B,C,D,E) = [ ,1, , ,1 , atu = -^. 

v ' ' ' ' ' V 100 100 100 10/ 100 



(4.23) 



There are two singularities located at Uh\ = 3.97 x 10 3 and Uhi = 1.27 x 10 l . Figure [3] shows that 
near Uhu D{u) and E(u) behave like 



and near Uhi, their behavior is 



D(u) ~ (u-Uhi) p , < p < 1 , 
E(u) ~ (u — Uhi) q ', 0<g<l, 



D(u) ~ (u h2 - u) p , p < 0, 
E(u) ~ (uhi — u) q , 0<q<l, 



(4.24) 



(4.25) 



with some positive numerical factors. Through an argument similar to that for the type I solution, we 
conclude that Uhi and Uhi are the event horizons. 
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Figure 4: Type III solution: plots of D and E. Horizontal lines are u. 



type III solution 

An example of the type III solution is given by the boundary condition: 

(A,B,C,D,E)= (— ,1-— ,— ,— ,1-— \ atu=4^- (4.26) 
V 100 10 100 100 10/ 100 v 

There is a singularity located at «/, = 4.1 1 x 10~ 3 . Near Uh\, D(u) does not have divergent behavior 
(at least numerically). It implies that D(u) and E(u) in figure [4]behave like 

D(u) C+D h (u-u h ) p , 0<p<l, 

E{u) ~ E h (u- u h ) 9 , 0<q< 1. (4.27) 

We conclude that F/G vanishes at «/, using an argument similar to the previous ones. Then u = u/, is 
the horizon. 

type IV solution 

An example of the type IV solution is given by the boundary condition: 

N , 1 1-2 1 1 . , 

(A,B,C,D,E) = ( ,1 , , ,H ), at« = -r^. (4.28) 

v 100 10 100 100 10 ; 100 

Numerical plots of this solution are displayed in figure [5] There is a singularity located at u < 10~ 20 . 
We think that the real singularity is located at u = 0, which is caused by the power law behavior of 
D and E near u = 0. There are no other singularities. Thus we interpret this type IV as the solution 
which describes a spacetime without black hole. 

Finally we point out that there are similarities between some of numerical solutions in this sub- 
section and the analytic solution (14. 10b - In fact the type III solution is similar to (14.101 ) with c < 0, 
and the type IV solution is similar to ( 14.101 ) with c > 0. 
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5 Summary and Discussion 

We have studied classical solutions of the torsion gravity in empty spacetime, whose field equation 
stems from the matrix equation of motion of the IKKT model by using a HKK interpretation. The 
PGT is one of the known torsion gravities which is formulated by a Lagrangian of the vielbein and 
the spin connection. Our field equation is not derived from such Lagrangian, so it is different from 
the PGT field equation, even though a matrix equation of motion itself is derived from a matrix model 
action. Because of the lack of the gravitational Lagrangian, no procedure for direct quantization of 
the vielbein and the spin connection is known in our approach. If their quantization rule is the same 
as the ordinary one, there is an anxiety about negative norm states caused from time-components 
of these fields as discussed in (H. In this paper we treat these gravitational fields and their field 
equations as purely classical objects to describe semi-classical aspects of the large N matrix model. 

Time dependent solutions with homogeneity and isotropy were investigated in section 3. We 
found that introducing torsion provides a simple mechanism to realize an accelerating expansion of 
spacetime. In PGT torsion degrees of freedom have been used to study cosmological problems like 
dark energy and inflation (see lfl4ll . for example). To carry out further investigation along this line, 
we need to take care of contributions corresponding to the energy-momentum tensor from matter, 
which is beyond our scope in this paper. 

Time independent solutions with spherical symmetry were investigated in section 4, with both 
analytical and numerical methods. We found solutions with singularities, which relate to the event 
horizons according to our criterion. Except for the formal power series solution, all solutions with 
the torsion are not asymptotically flat. The asymptotic behavior is power law. It is contrast to PGT, in 
which an asymptotically flat solution has been found analytically lfT5l (several asymptotically non- 
flat solutions were also found iflolO . This would imply that torsions in our study are constrained to 
be very small in physically acceptable solutions, in order to well describe the solar system in present 
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dajO However this observation itself does not prohibit a possibility of finding some solutions with 
a large effect of torsion in past or future period, which describes possible evolutions of the universe 
. More precise investigation on the structure of these solution and their systematic classification are 
future problems. 
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A Bianchi identities 

We use differential forms T c = \T ab c Q a A B b and R c d = \R a b c d^ A b with 6" = e^chtf. Then (l24l) 
can be written as the following form: 

T a = dQ a -(i) b a AQ b , (A.l) 
R a b = d(o"b + (o" c A (Ofo c (A.2) 

where exterior derivative d acts as dF = dx^e^'d^ AF on a differential form F . We have introduced 
a one form (£> a b = G) CM b Q c . To obtain CO, we also used a property e^d a e^ = —e^d^ which follows 
from e c /J e' J b = § c z>- 

Now we may drive Bianchi identities using (IA. lb and (IA.2b . First we take the exterior derivative 
of dA.ll ), then obtain 

= dT d + dm x d AQ X -(o x d AdQ x 

= 2 ( d xT YZ d e x a e y a e z + r l YZ dQ Y a e z - t cI yz e F a je z ) + dxo x d a e x - ® x d a d$ x 

(A3) 

Using dA.ll ) and dA.21 ), then dA.31 ) becomes 

= ^d x T d YZ Q X f\Q ¥ f\$ Z + ^T d YZ ((O s Y Q s AQ Z -Q ¥ A (O s z Q s )-(O s d AT S 
+ br d YZ {T Y Ae z -e 5 ' AT Z ) 

+/? s d a e 5 

-(Ox 5 A (0 s d A e x - ©x^ A ©/ A e 5 

= ^(3 fl r**c + r rf yc (o fli / - r rf uZ (o fo , c z - (o a /T s bc + r^r 1 ^ +/w , )e fl a e b a e c 

= i iy a T d bc + T d YcT Y ab + RalJW A 8* A V . 
9 Actually torsions in the standard model extension are constrained to be order 10~ 31 Gev by measurement |17| . 
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In this way we obtained the first Bianchi identity 

V(aT d , h c) + T d Ac T e >ab) +R {abc) d = 0. (A.4) 
Similarly, from (1A.2I) . we can derive the second Bianchi identity 

V ia R bc) d e -Tf {tab R c) / e = 0. (A.5) 



Here we see that equations ( 12.91 ) which follow from the matrix Jacobi identity are identical with 
Bianchi identities (|A~4b and (|A3T ). Thus all T c \ ab and R ab cd defined by (H3]> satisfy (|2T9l . 
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